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I  INTRODUCTION 


Laser-induced  fluorescence  (LIF)  provides  a  sensitive  means  of  detecting 
trace  species,  often  free  radicals,  that  are  the  reactive  intermediates  in 
combustion  chemistry.*  Qualitative  identification  can  provide  insights  into 
the  mechanism  of  that  chemistry,  while  quantitative  measurement  yields  data 
for  comparison  with  computer  models  containing  detailed  kinetics.  One  of  the 
attributes  of  LIF  is  its  ability  to  provide  spatially  and  temporally  resolved 
information.  Point-resolved  measurements  are  made  by  focusing  the  fluores¬ 
cence  at  right  angles  to  the  laser  beam  onto  a  slit  oriented  perpendicular  to 
the  beam  propagation  direction,  whereas  temporal  resolution  is  Inherent  in  the 
use  of  pulsed  lasers.  This  has  been  the  customary  mode  of  acquisition  of  LIF 
data  on  stable  flames  where  the  burner  is  translated  through  the  optical  mea¬ 
surement  point,  so  as  to  yield  a  profile  of  the  pertinent  species. 

In  a  time-varying  flow  such  as  turbulent  combustion,  however,  the  system 
changes  character  at  a  given  point  between  successive  laser  shots  (which  are 
typically  separated  by  100  msec  for  lasers  used  in  such  studies.)  Thus  a 
series  of  point  measurements  can  be  treated  only  statistically,  in  the  form  of 
probability  distribution  functions,  for  comparison  with  theoretical  treat¬ 
ments.  It  would  instead  be  highly  desirable  to  obtain  correlated  information 
during  a  single  laser  pulse.  Such  Information  could  take  two  forms:  The 
spatial,  multipoint  distribution  of  some  observable,  or  the  measurement  of  two 
or  more  observables  simultaneously.  The  two  could,  in  principle,  then  be  com¬ 
bined  to  determine  the  instantaneous  spatial  distribution  of  a  pair  of  observ¬ 
ables. 

The  demonstration,  development,  and  study  of  methods  to  meet  these  needs 
has  been  the  objective  of  this  project.  The  instanteous  spatial  distribution 
of  OH  molecules  in  a  flame  can  now  be  determined  using  the  technique  of 
two-dimensional  fluorescence  imaging,  first  performed  by  us  in  this  work^*^ 
and  now  in  use  in  several  laboratories  around  the  world.  The  simultaneous 
measurement  of  OH  concentration  and  the  temperature  has  been  addressed  using 
the  method  of  vibrational  energy  transfer,^  although  we  have  not  been  success- 


ful  in  determining  a  foundation  for  implementing  the  method  quantitatively  for 
a  range  of  flame  conditions.  The  OH  radical  has  been  the  subject  of  these 
studies.  Ubiquitous  in  flames  and  important  in  many  of  the  major  and  minor 
chemical  pathways,  the  OH  radical  is  often  present  in  relatively  large  amounts 
(>  100  ppm)  and  its  absorption  lies  in  an  experimentally  convenient  laser 
wavelength  region.  It  is  the  key  intermediate  detectable  by  laser  methods  and 
the  one  most  likely  to  be  of  Interest  in  measurements  on  turbulent  flames. 

II  TWO-DIMENSIONAL  FLUORESCENCE  IMAGING 

Our  efforts  on  the  project  were  devoted  to  the  demonstration  and  develop¬ 
ment  of  two-dimensional  fluorescence  Imaging.  This  proceeded  quickly  and 
successfully,  yielding  the  first  signals  in  January  1982.  Subsequent  study 
and  refinement  of  the  method  showed  excellent  sensitivity  on  single  laser 
shots.  Following  our  announcement  of  the  demonstration,  the  method  was 
successfully  attempted  elsewhere  and  is  now  in  use  in  at  least  a  half-dozen 
laboratories.  The  method  and  the  results  of  our  study  of  it  were  fully 
described  in  two  paper s^»^. 

Ill  VIBRATIONAL  ENERGY  TRANSFER  THERMOMETRY  (VETT) 

The  objective  of  this  phase  of  the  work  was  the  development  of  a  method 
to  provide,  on  a  single  laser  shot  with  one  laser,  a  simultaneous  measurement 
of  the  OH  concentration  and  the  temperature.  The  concept  of  the  method^  is  as 
follows. 

Figure  1  depicts  the  energy  levels  and  the  radiative  transitions  involved 
in  the  experiment.  A  laser  is  tuned  to  the  wavelength  of  a  given  absorption 
line  of  the  (0,0)  band  of  the  A^E+-X^n^  electronic  transition,  populating  a 
particular  rotational  level  in  v'  ■  0.  (The  rotational  levels  are  not  shown 
explicitly  in  Figure  1.)  While  in  the  A^E  +  state  and  before  radiating,  a  few 
of  the  molecules  excited  by  the  laser  undergo  collisions,  raising  them  into 
the  v'  ■  1  vibrational  level.  Measurements  of  the  intensities  of  the  (0,0) 
and  (1,0)  bands,  as  indicated  in  the  figure,  furnish  the  ratio  of  the  popula- 


tions  in  the  two  vibrational  levels,  from  which  the  temperature  can  be 
derived. 


The  population  ratios  reported  in  Ref.  4  were  analyzed  within  the  frame¬ 
work  of  a  steady-state  approximation.  We  consider  only  the  total  population 
in  each  vibrational  level;  that  is,  any  influence  of  rotation  on  the  vibra¬ 
tional  transfer  collisional  processes  is  ignored.  A  steady-state  balance 
applied  to  Nj  yields 


[V  exp(-AE/kT)]NQ  -  (V  +  A  +  Q)Nt  .  (1) 


Here  (see  Fig.  3),  the  downward  vibrational  transfer  rate  is  V,  and  the  upward 
transfer  is  given  by  detailed  balancing  as  V  exp(-*\E/kT).  Q  is  the  quench 
rate,  and  A  is  the  Einstein  emission  coefficient  for  v’  ■  1.  Since 
A/(V  +  Q)  <  0.01,  A  was  ignored.  Rearranging, 


N1  _  exp(-AE/kT) 

N  ’  1  +  Q/V 

o 


(2) 


.  N 

T - £-  In  (1  +  Q/V)  -2-  .  (3) 

o 

Thus  from  the  measured  ratio  and  a  value  of  Q/V,  one  may  determine  the 
temperature. 

The  method  was  originally  demonstrated^  using  point  measurements  in  the 
burnt  gases  of  (Hfy/air  flames.  Temperatures  Ty  measured  by  VETT  were  found  to 
agree  well  with  measurements  of  the  ground  state  rotational  temperatures  TR 
obtained  via  excitation  scans  (after  subsequent  corrections^  to  each  were 
applied,  necessary  because  of  the  neglect  of  rotational  transfer  effects**  on 
Tr  and  an  incorrect  estimate^  of  Q/V  in  Ref.  4).  Good  agreement  (within  10%) 
was  achieved  between  Ty  measured  following  excitation  of  three  different 
rotational  levels  in  v'  ■  0  (F1(4),  F2(10),  F2 C14) ] .  Using  split  images  and 
separately  filtered  vidicon  halves,  we  could  in  principle  use  this  method  in 
an  imaging  experiment  to  yield  an  instantaneous  map  of  the  OH  concentration 
and  temperature  throughout  those  regions  of  the  flame  where  OH  was  present. 


Following  the  demonstration  of  VETT  in  Ref.  4,  several  questions  remained 
concerning  its  applicability  in  a  quantitative  way,  for  a  range  of  tempera¬ 
tures  and  flame  gas  collisional  environments.  These  are 

(1)  Given  the  knowledge  of  rotational-state-specific  vibrational  energy 
transfer  in  A^£+OH,  both  in  room  temperature  flows^  and  in  flames, 
what  are  the  effects  of  pumping  different  N’  levels  in  v'  «  0  and 
measuring  only  part  of  the  rotational  distribution  in  v'  ■  1? 

(2)  How  does  the  parameter  Q/V  (needed  for  the  analysis)  vary  with  the 
collisional  environment,  both  gas  mixture  and  temperature? 

(3)  Over  what  range  of  temperature  can  one  reliably  measure  Ty  and  what 
1 8  its  precision,  compared  to  nonintrusive  laser  measurements  of  TR 
as  well  as  thermocouple  measurements  that  furnish  Ttc  but  are 
intrusive  and  require  radiative  corrections? 


These  questions  have  been  addressed  by  a  series  of  single-point , 
multi-pulse  averaged  measurements  in  the  burnt  gases  of  stable  flames  of 
CH^/02/air.  These  flames  have  furnished  a  temperature  range  of  about  1400  to 
2700K  with  high  enough  OH  densities  to  produce  usable  signals. 


A.  Signal  Analysis  Equations;  Determination  of  Q/V 

Eq .  (3)  was  derived  under  two  assumptions:  first,  a  steady  state  treat¬ 
ment  of  the  populations  and,  second,  that  the  quenching  rate  Q0  from  v'  ■  0  is 
the  same  as  the  quench  rate  Qj  from  v'  ■  1  (each  denoted  by  Q).  We  here 
consider  the  time-dependent  equations  that  describe  the  actual  experiment. 

Neglecting  rotational  level  effects,  the  pertinent  equations  are:  (i) 
Uptransfer,  pumping  v'  *  0: 

dN 

d r  “  L(t)  ‘  QoNo 

Ve^E/kT  N0  -  QjNi  -  VHj  (5) 


R 


u 


rQ  Nl(t)dt/r0  N0(t>dt 


(6) 


(ii)  Downtransf er,  pumping  v'  ■  1 : 


d r  “  L(t)  "  VNi  “  QiNi 


dT  -  VN1  -  ^oNo 
RD  -  C  N0(t)dt/r  Nj(t)dt 


Here,  L(t)  is  the  time-dependence  of  the  laser  pulse;  it  is  conveniently  and 
realistically  exjr^ssed  as  the  difference  of  two  exponentials: 

L(t)  ■  e  *  -e  ^  .  In  the  experiment,  we  use  a  gated  boxcar  integrator  whose 
output  is  proportional  to  the  intensity  of  some  band  [(1,0),  (1,1)  or  (0,0)] 
from  which  is  obtained  the  population  of  the  corresponding  v*  level.  The 
detector  time  constants  are  much  longer  than  the  inverses  of  Ip  or  the 
collisional  rates  involved.  Thus  the  boxcar  output  yields  f*  N^(t)dt  and  the 
measured  band  ratios  are  those  given  by  Ru  and  Rp. 

An  equation  of  general  form 


\  exp  (-ait)-  cx 


is  solved  by 


l 

1  _  — .  (  exp(-a  t)  -  exp(-ct)J 

i  C  a  .  1 


and  has 


r  xdt  -  -  e  — 

o  c  J  at 


Therefore, 


'D  V 


V  exp  (-AE/kT) 
^  *  V  +  Ql 
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Note  that  Q  the  quenching  rate  for  v'  ■  0,  is  determined  by  the  downward 
transfer  experiment  pumping  v*  ■*  1,  whereas  Q^,  the  quenching  rate  for  v'  *  1, 
is  needed  for  analysis  of  the  v'  -  0  pump  uptransfer  runs  used  to  determined 
the  temperature. 


Experimental  Arrangement 


The  experiments  described  in  the  succeeding  sections  were  all  performed 
using  point-resolved  measurements  in  the  burnt  gases  of  premixed  CH4/O2/N2 
flames  at  atmospheric  pressure.  The  burner  was  1  cm  in  diameter;  its  surface 
consisted  of  a  plate  drilled  with  1-mm-d  holes  to  ensure  even  flows  in  the 
burnt  gas  region.  The  horizontal  laser  beam,  1  mm  in  diameter,  passed  across 
the  burner  and  the  fluorescence  was  focused  at  right  angles  onto  the  vertical 
slit  of  a  3/4-tn  spectrometer.  Slit  widths  ranging  from  0.1  to  3  mm  were  used, 
depending  on  the  experiment  conducted.  The  fluorescence  was  detected  with  an 
EMI  9558QA  photonultlplier.  In  a  few  excitation  scan  runs,  total  undispersed 
fluorescence  was  collected  with  a  filtered  1P28A  photomultiplier  viewing  the 
center  of  the  flame.  The  photomultiplier  output  in  each  case  was  fed  to  a 
boxcar  Integrator.  The  boxcar  output  in  turn  was  recorded  on  a  strip  chart 
recorder  for  some  runs;  for  other  runs  the  output  was  stored  in  the  Molecular 

Q 

Physics  Laboratory  computer  through  a  remote  hookup.  The  flames  were 

generally  slightly  rich  ($=1.1  to  1.2)  mixtures  of  CH^ ,  C>2 ,  and  air.  The 

flows  were  measured  with  rotameters  that  had  been  calibrated  with  a  wet  test 

q 

meter.  According  to  equilibrium  calculations,  the  OH  concentration  varied 
from~  30  to  ~  1500  ppm  over  the  series  of  mixtures  used;  occasionally, 
absorption  measurements  were  made  that  confirmed  these  values.  Tests  were 
made  to  ensure  operation  in  regions  of  the  flame  where  absorption  by  OH  of  the 
laser  light  or  the  fluorescence  did  not  pose  problems. 


C.  Rotational  Distributions  for  Uptransferred  OH 

In  this  phase  of  the  study,  we  measured  the  rotational  distributions  of 
OH  molecules  that  had  been  collisionally  transferred  upwards  from  v'  ■  0  to 
v'  ■  1.  Individual  rotational  levels  within  v'  =  0  were  pumped  by  the  laser, 
and  rotationally  resolved  fluorescence  scans  were  made  of  the  (1,0)  band  in 
fluorescence  to  determine  relative  populations  among  the  rotational  levels  in 
v'  “  1.  This  information  is  needed  to  ascertain  the  spread  of  the  spectrum. 


6 


r 


to  set  wavelength  observation  limits  for  the  less  well-resolved  thermometry 
runs,  to  choose  the  pertinent  N'  In  v'  ■  1  that  should  be  pumped  to  measure  a 
pertinent  Q/V  In  the  downward  transfer  measurements,  and  to  address  questions 
concerning  possible  N'  variation  of  both  Q  and  V  that  may  enter  Into  the  data 
Interpretation. 

Using  a  variety  of  slit  widths  furnishing  resolution  between  0.2  and  1  A  , 
we  obtained  fluorescence  spectra  In  the  burnt  gas  regions  where  excitation 
scans  showed  the  rotational  temperature  TR  ~  2100K.  A  sample  spectrum  Is 
shown  in  Figure  2.  From  such  spectra,  50  spectral  "features”,  definable  peaks 
or  shoulders,  were  Identified  and  their  intensities  measured.  Each  of  these 
features  is  composed  of  from  1  to  8  known  rotational  lines  of  the  (1,0)  band, 
the  exact  mix  depending  on  the  spectrometer  bandwidth  used  for  the  particular 
run.  The  intensities  are  then  Interpreted  in  terms  of  relative  populations  of 
the  31  rotational  levels,  F^N'  )  and  F2(N' )  for  N'  ■  0  to  15.  An  analysis 
program  was  written  using  listed  line  positions10  and  line  strengths.11  The 
spectrometer  bandwidth  AX.  could  be  varied  over  appropriate  limits,  and  a  least 
squares  fit  to  the  populations  was  made.  For  example,  in  some  of  the  earlier 
lower  resolution  runs,  a  variation  of  the  assumed  AX  between  0.80  and  0.95  A 
reduced  the  X  value  by  30%.  Most  of  the  later  runs  were  made  with  higher 
resolution,  near  0.2  Af  where  only  1  or  2  levels  contributed  significantly  to 
each  chosen  feature. 

The  results  from  pumping  a  series  of  levels  are  shown  in  Figure  3  in  the 
form  of  Boltzmann  plots  for  ready  comparison  with  TR.  There  is  some  scatter 
at  the  lowest  N' ,  due  to  low  intensities  for  the  pertinent  features  and  conse¬ 
quent  ill-conditioning  In  the  fit;  this  scatter  does  not,  however,  affect  the 
essential  conclusions  and  no  real  effort  was  made  to  improve  the  fitting 
procedure  to  accommodate  the  problem. 

We  may  conclude  that,  within  experimental  error,  the  population  distribu¬ 
tion  of  the  uptransf erred  molecules  is  near-Boltzmann  in  nature,  described 
reasonably  by  a  temperature  very  close  to  the  gas  temperature,  and  has  no 
significant  dependence  on  the  originally  pumped  level  in  v'  -  0.  Note  that 
the  N'  *  5  and  6  levels  have  but  500-700  cm-1  of  rotational  energy,  whereas 
N'  ■  12  and  13  are  at  2600  and  3000  cm-1  above  N'  ■  0  of  v'  ■  0;  these  latter 
energy  levels  may  be  compared  with  AGj/2  of  the  A  £  state,  which  is 
2990  cm"1. 
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These  results  in  turn  suggest  that  dovntransfer  experiments  pumping 
H'  ■  6  in  v'  «  1,  which  is  the  most  highly  populated  rotational  level  in  this 
temperature  region,  should  well  simulate  the  collislonal  behavior  of  the  ufr- 
ward  transferred  OH,  assuming  that  enough  rotational  redistribution  takes 
place  following  the  Initial  excitation  of  N'  ■  6. 

D.  Upward  Transfer  as  a  Function  of  Rotational  Level  in  v'  ■  0 

In  the  initial  VETT  experiments,^  measurements  of  the  degree  of  upward 
0  *  1  transfer  were  made  by  pumping  different  rotational  levels.  The  results 
for  N'  ■  4,  10,  and  13  showed  the  same  degree  of  transfer  within  experimental 
error.  This  was  surprising,  for  it  might  have  been  anticipated  that  the 
additional  rotational  energy  would  facilitate  the  endoergic  transfer.  (Down¬ 
ward  transfer  has  been  shown'’  to  decrease  with  increasing  H'  in  v'  *  1).  Ue 
have  carried  out  a  more  definitive  set  of  scans,  in  two  different  flames,  and 
find  results  conflicting  with  those  of  Ref.  4.  Here,  the  upward  transfer 
ratio  of  Nj/Nq  was  measured  and  found  to  vary  smoothly  with  N' .  We  do  not 
understand  the  disagreement  with  the  results  of  Ref.  4,  although  a  slightly 
different  flame  environment  was  used  in  the  current  scans. 

The  first  set  of  runs  was  made  on  downward  transfer  to  check  the  correc¬ 
tion  for  band  overlap.  The  (1,1)  band  head  is  at  3122  A.  p  branches  of  the 
(0,0)  band  for  N'  >  8  and  Q  branches  for  N'  >  15  are  to  the  red  of  this  wave¬ 
length.  Thus  the  measured  intensity  between  3064  and  3122  A  represents  only  a 
fraction  f  of  the  true  (0,0)  intensity,  while  that  from  3122  A  onward  to  the 
red  Includes  all  of  the  (1,1)  band  plus  1-f  of  the  (0,0)  band.  Since  the 
fraction  of  molecules  with  quantum  number  JD  or  greater  at  some  temperature  T 
is 


J“  (2J  +  l)e-BJ(J  +  1)/kTdJ/f  (2J  +  l)e“BJ(J  +  1)/kTdJ 


(15) 


-BJ  (J  +  l)/kT  , 
o  o 


f  may  be  evaluated;  for  a  2200K  thermal  distribution  it  is  0.85.  Using  this 
correction  for  the  (1,1)  band,  measurements  were  made  of  the  v'  *  1  population 
using  both  the  (1,0)  and  the  (1,1)  band  intensities.  Needed  are  the  spectro- 


meter  and  photomultiplier  responses  at  the  respective  wavelengths;  these  are 
taken  from  manufacturers'  curves  for  our  EMI  tube  and  Bausch  and  Lomb 
grating.  Also  needed  are  the  relative  Einstein  coefficients  for  the  two 
bands. The  results  for  the  two  bands  yield  N^  "  22.08  and  21.95,  agreeing 
better  than  the  estimated  intensity  uncertainties  and  indicating  that  the 
correction  is  properly  accounted  for. 

Measurements  of  the  upward  transfer  ratio  N^/Nq  were  then  made  in  two 
separate  flames:  cool  (1650K)  and  hot  (2800K),  where  the  temperatures  were 
determined  by  excitation  scans.  Exemplary  data  are  shown  in  Figure  4.  In 
each  case  three  different  rotational  levels  were  pumped.  The  results  are 
given  in  Table  1.  Note  that  the  two  Qj9  runs  for  the  cool  flame  and  the  Q^9, 
PjlO  pair  for  the  hot  flame  constitute  replicate  runs  because  in  each  case  the 
same  upper  level,  F^(9),  is  excited.  The  agreement  for  these  pairs  is 
excellent. 

In  each  case  the  hotter  flame  shows  more  upward  transfer,  as  expected. 
However,  the  Nj/Nq  ratio  increases  smoothly  with  increasing  N'  in  v'  ■>  0,  in 
contrast  to  the  results  of  Ref.  4.  We  did  not  attempt  to  address  this  discre¬ 
pancy.  Neither  did  we  pursue  the  interesting  fundamental  questions  posed  by 
these  results  concerning  the  efficacy  of  rotational  energy  in  promoting  upward 
vibrational  transfer  at  the  expense  of  large  changes  in  angular  momentum. 
Rather,  we  decided  to  conduct  the  remainder  of  the  experiments  pumping  N*  ■  6, 
J*  -  61/2,  In  both  v'  ■  0  for  upward  transfer  and  v'  ■  1  for  downward 
transfer.  This  is  the  most  highly  populated  level  for  thermal  populations 
near  2000K  and  in  the  rotational  distribution  of  upward  transferred  molecules, 
as  described  in  the  preceding  section.  The  behavior  of  the  N'  -  6  level 
should  therefore  be  reasonably  considered  as  that  of  an  average  N'  value, 
given  a  smooth  variation  in  collision  rate  with  N' ,  and  the  effects  of  any 
rotational  redistribution  that  does  occur  will  be  minimized. 

E .  Final  Q/V  and  Thermometry  Runs 

In  this  section  we  describe  our  final  set  of  runs  measuring  T,  Q/V,  and 
Ni/N0  in  a  series  of  flames. 


(1)  Flames  Used.  A  series  of  three  flames  were  probed  so  as  to  produce 
different  but  overlapping  temperature  ranges,  although  each  had  a  slightly 


‘'j'* ^  .*■  .  *  A 
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Table  1 


N1/Nq  RESULTS  FOR  UPWARD  TRANSFER  FOR  DIFFERENT  PUMP  LINES 


LI  ne 

Cool  flame 

Hot  flame 

P21* 

0.0327 

0.0417 

V 

0.0406 

- 

Qi9 

0.0456 

0.0498 

P^O 

- 

0.0499 

Qx13 

0.0527 

0.0612 

1  0 


varying  collislonal  mix  in  the  burnt  gas  region.  The  flows  of  methane, 
oxygen,  and  air  used  are  given  in  Table  2;  we  shall  denote  the  flames  as  H,  M 
and  L  in  the  following,  in  spite  of  the  temperature  overlap.  Operation  at 
different  points  (heights  above  the  burner)  in  the  burnt  gases  yielded  a  wide 
range  of  usable  temperatures,  from  1400  to  2700K. 

(2)  Temperature  Determination  by  Thermocouples  and  Excitation  Scans*^. 

We  devoted  considerable  effort  to  ensuring  that  we  could  accurately  determine 
the  gas  temperature  at  different  points.  First,  a  series  of  excitation  scans 
was  made  to  determine  TR.  Early  replicate  runs  showed  disagreement  outside 
fitting  errors,  but  this  was  found  to  be  due  to  optical  depth  problems.  A 
typical  Boltzmann  plot  from  one  of  the  later  runs  is  shown  in  Figure  5. 

Except  for  the  lowest  P-branches,  the  data  fit  with  standard  deviations  of 
3-5Z.  The  P-branch  discrepancy  may  be  due  to  polarization  effects*^  but  the 
cause  was  not  investigated  here,  and  the  Q-branch  sequence  was  used  to 
determine  TR. 

An  extensive  series  of  thermocouple  measurements  was  made  in  each 
flame.  Four  thermocouples  of  varying  bead  size  (nominal  diameter  0.002", 
0.005",  0.015",  and  0.020")  were  used  to  provide  related  values  of  the  radia¬ 
tive  correction.  Measurements  were  made  at  as  many  as  40  points  in  each  flame 
for  each  thermocouple,  at  different  positions  through  the  burnt  gas  region  as 
close  to  the  flame  zone  as  survival  of  the  thermocouple  would  permit.  A 
series  of  measurements  is  shown  in  Figure  6. 

A  radiative  correction  is  necessary  because  of  emissive  heat  loss  from 
the  thermocouple.  The  gas  temperature  Tg  is  related  to  the  bead  temperature 
by 

oeD  (T.4  -  T  4) 

Ts  '  Ts  + - kW  —  «‘> 


where  a  is  the  Stefan-Boltzmann  constant,  D  and  t  the  bead  diameter  and 
emissivlty,  k  the  gas  thermal  conductivity,  Nu  the  Nusselt  number,  and  Ty  the 
temperature  of  the  radiation  sink  (300K).  Using  the  nominal  bead  diameters 
and  c  *  0.3,  we  were  unable  to  obtain  consistency  among  Tg  values  calculated 
for  each  thermocouple  at  each  point.  Therefore,  Eq.  (16)  was  used  in 
empirical  form  (noting  Tw^  <<  T^): 


Table  2 


FLOW  RATES  (cm3  sec-1)  AND  OTHER  PARAMETERS  FOR  FLAMES  IN  FINAL  RUNS 


H 

M 

L 

ch4 

12.7 

9.78 

6.52 

°2 

13.5 

11.9 

0 

Air 

42 

29.2 

52.5 

4> 

1.14 

1.09 

1.18 

Temperature  Range 

1900-2700 

1600-2400 

1300-2000 

Used  (K) 


(2) 


A  single  value  of  a  was  determined  for  all  the  measurements  such  that 
differences  in  Tg  for  the  0.002",  0.005",  and  (in  the  H  and  M  flames)  0.015” 
thermocouples  were  minimized  at  each  point.  (The  0.020”  set  was  not  consis¬ 
tent.)  The  values  of  a  so  determined  were  aR  *  1.33  x  10~8,  Sju  -  1.25  x  10~8, 
a^  •  1.51  x  10~8  compared  to  a  calculated  value  0.66  x  10-8  from  the  constants 
of  Eq.  (16),  all  in  units  of  cm~^  K-^.  The  difference  reflects  primarily  a 
larger  effective  bead  diameter  compared  to  the  nominal  value.  A  value 

_  Q 

a  *  1.32  x  10  was  used  to  reduce  the  remainder  of  the  data.  The  resulting 
temperatures  then  agreed  to  within  1-3 X.  Results  for  the  VETT  measurement 
locations  are  given  in  Table  3.  The  highest-T  values  from  the  0.015" 
thermocouple  were  not  used  when  they  clearly  disagreed  with  those  of  the 
smaller  bead  devices.  Values  in  parentheses  are  extrapolations  from  lower 
temperature  along  curves  parallel  to  those  of  the  larger  thermocouples. 

The  results  also  agreed  with  excitation  scan  values  of  TR  at  the  same 
point.  For  example,  at  the  3.1  and  3.6  cm  positions  of  flame  H,  we  determined 
Ttc  -  2068  ±  34K,  Tr  -  2032K,  and  Ttc  -  2141  ±  38K,  TR  -  2157K, 
respectively.  The  uncertainties  in  Ttc  represent  the  spread  in  values 
obtained  by  applying  Eq .  (17)  to  the  three  thermocouples,  all  with 
a  -  1.32  x  10*8. 

In  the  following  sections,  Q/V  and  Nj/N0  measurements  are  described. 

These  were  made  at  a  total  of  20  spatial  locations  (6  in  flame  L,  9  in  M,  5  in 
H).  In  each  case  N'  *  6,  J  ■  6.5  level  was  pumped  in  both  v'  •  0  and  1.  The 
highest  temperature  point  in  L  was  clearly  erroneous,  because  the  finite-size 
laser  beam  partially  covered  a  region  of  sharp  temperature  gradient  in  the 
reaction  zone,  and  it  is  excluded  from  further  discussion.  The  other  19 
measurements  are  all  included.  At  each  spatial  location,  Ttc  was  determined 
using  the  procedure  described  in  this  section.  This  Ttc  forms  the  value 
plotted  on  the  abscissa  of  Figures  7-10. 


(3)  Q/V  Results.  The  value  of  Qq/V  was  determined  from  the  degree  of 
downward  transfer  (see  Eq .  13).  The  method  of  Ref.  5  was  used  to  analyze  the 
population  ratios  from  the  intensities  of  the  (1,0)  and  (0,0)  bands, 
correcting  the  latter  for  the  overlap  with  (1,1)  as  described  above.  The 
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Table  3 

THERMOCOUPLE  RESULTS  (K)  FOR  a  -  1 .32  x  10-9 
(Asterisk  denotes  not  used;  parentheses  denote  extrapolated  value) 

_ at _  _ t _ 

Flame  and  Point  0.002”  0.005"  0.015"  0.002"  0.005"  0.015 


results  are  listed  In  Table  4  and  plotted  In  Figure  7  as  a  function  of  Tt  for 
the  three  flames. 

A  very  large  variation  was  seen  because  the  v*  *  0  rotational  distribu¬ 
tion  may  vary  with  gas  temperature,  and  collislonal  quenching  rates  from 
v*  ■  0  are  known  to  vary  with  N* ,  according  to  measurements  at  room 
temperature. 1^-17  Additionally  the  collislonal  environment  varies  slightly 
from  flame  to  flame,  although  not  enough  to  change  Qq  as  much  as  is  seen  here, 
according  to  measurements  at  room  temperature1 ^  and  elevated  temperature1®  for 
a  variety  of  collision  partners.  V  also  varies  with  collider  gas  according  to 
room  temperature  data,^*1^  although  vibrational  transfer  of  A%+  OH  has  never 
been  measured  for  ^0  as  a  collision  partner.  The  fact  that  all  the  points 
for  a  given  flame,  H,  M  and  L,  lie  on  a  smooth  line  —  which  is  different  for 
each  flame  —  suggests  that  differences  in  the  collislonal  mixture  is  the 
cause  of  the  Q/V  variation.  Nonetheless,  the  variation  of  Q/V  is  much  greater 
than  expected  on  these  grounds. 

W  VWq  results .  Upward  transfer  runs  were  made  at  the  same  spatial 
locations  of  the  Q/V  measurements.  Analysis  was  as  described  in  Ref.  4  and 
Section  IIIA;  the  results  are  listed  in  Table  4  and  are  plotted  versus  Ttc  in 
Figure  8.  Note  that  there  is  a  smooth  overall  increase  in  N1/NQ  with 
increasing  temperature,  as  anticipated. 

F.  Thermometry  Results 

(1)  VETT  Results.  Eq .  (3)  was  used  to  determine  the  temperature  Tv  at 
each  measurement  point,  taking  as  input  the  experimental  values  of  RD  *  Q/V 
and  Nj/Nq.  The  results  are  listed  in  Table  5  and  plotted  versus  Ttc  in 
Figure  9.  They  do  not  make  sense.  Although  N^/Nq  increases  with  temperature 
as  expected,  the  increase  in  Q/V  with  Tt  (Figure  7)  forces  Tv  to  rise  much 
too  sharply  and  attain  wholly  unrealistic  values. 

(2)  An  Approach  Using  an  Effective  Value  of  Q/V).  The  lack  of  correspon 
dence  of  Ty  with  Tfcc  may  be  due  to  rotational  dependences  in  the  collision 
rates,  or  the  fact  that  Q0/V  is  measured  by  the  downward  transfer  (Eq .  13), 
whereas  Q^/V  is  needed  to  obtain  Tv  (Eq.  14).  The  fact  that  relatively  smooth 
variation  is  seen  for  each  flame  (Figures  7  and  8)  encouraged  us  to  try  a 
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Table  4 

RESULTS  OF  THERMOMETRY  MEASUREMENTS 


Measurement 


Flame 

Point 

Ttc 

Nj/N0x  102 

Q/V 

Tv  (Eq.  3) 

L 

1 

1809 

2.27 

1.83 

1519 

2 

i.712 

2.26 

1.70 

1491 

3 

1589 

2.14 

1.60 

1443 

and  the 


different  approach.  The  temperature  at  each  point  was  equal  to  Tt  , 
measured  ratio  was  then  used  to  derive  an  effective  Q/V  by  inverting  the 
thermometry  equation  (3): 


(-AE/kT  ) 
expw  t  c 

X7n 

1  o 


(18) 


The  results  from  this  step  are  given  in  Table  5. 

Table  5  shows  a  definite  ordering  of  the  values  of  (Q/V)eff  by  flame.  We 
thus  considered  that  this  parameter  may  be  the  same  for  all  the  points  in  a 
given  flame,  even  though  it  may  differ  from  flame  to  flame  because  of  a 
different  mixture  of  collision  partners.  Accordingly,  average 
values  (Q?V)  were  calculated.  (Only  6  of  the  9  points  in  the  M  flame  were 
used;  those  not  included  in  computing  the  average  are  marked  with  an  asterisk 
in  Table  5.)  The  results  are  (Q'/v)^  -  2.57  ±  0  62  for  L,  1.88  ±  0.12  for  M, 

and  1.47  ±  0.07  for  H.  These  are  far  from  the  values  measured  by  downward 
transfer  and  are  in  fact  in  the  opposite  direction  (compare  Tables  4  and  5). 

We  do  not  know  the  reasons  for  these  differences. 

Nonetheless  we  proceeded,  applying  the  average  (Q/v)  for  each  flame  to 

ef  f 

the  measurements  at  all  the  points  in  that  flame.  With  this  a  new, 
"calibrated”  Ty  (which  we  label  Tc)  was  calculated  using  Eq.  (3).  The  results 
are  listed  in  Table  5  and  plotted  versus  Ttc  in  Figure  10.  Note  that  the 
scale  is  different  from  that  in  Figure  9. 

These  "calibrated"  thermometry  results  are  in  very  good  agreement  with 
Ttc,  as  seen  by  the  straight  line  in  Figure  10  which  is  the  equation 
Tc  “  Ttc.  For  the  19  points,  the  standard  deviation  is  3.56%  and  rms 
deviation  is  4.96%.  In  fact,  4  points  with  a  deviation  of  8%  or  more  in 
(Tc-Ttc)/Ttc  contribute  about  half  of  this  uncertainty.  In  situ, 
spectroscopic  flame  temperature  measurements  with  this  level  of  accuracy  can 
be  considered  excellent. 


G.  Summary  and  Discussion  of  Thermometry  Results 

In  this  study  we  first  determined  that  we  can  make  precise  measurements 
of  the  temperature  in  the  burnt  gases  of  a  methane/air  flame,  using  rotational 
excitation  scans  and  thermocouple  measurements.  This  turned  out  to  be  a  more 


Measurement 

"CALIBRATED" 

Table  5 

THERMOMETRY  RESULTS 

Point 

T 

‘tc 

(Q/V)  Tv  (calibrated) 

Points  not  Included  In  computing  average  values  of  (Q/V)eff  for  the  M  flame 


painstaking  task  than  first  assumed.  There  existed  ample  opportunity  for 
systematic  errors  in  the  TR  determinations,  and  the  deviation  of  the  low 
P-branches  from  the  Boltzmann  plot  (Figure  5)  indicates  that  we  have  not 
solved  all  of  them.  Only  by  fitting  a  series  of  thermocouple  results,  using 
varying  bead  sizes,  were  we  able  to  attain  a  self-consistent  set  of  radiative 
corrections.  We  are  wary  of  any  flame  thermocouple  measurements  made  with  a 
single  thermocouple  using  the  nominal  equation  and  constants  of  Eq .  (16). 

The  VETT  temperature  Ty  was  then  determined  using  both  upward  and 
downward  transfer  to  determine,  at  the  same  locations,  N^/NQ  and  Q/V. 

Combining  these  to  obtain  Ty  produced  unrealistic  results  compared  with  the 
thermocouple  measurements.  The  results  make  too  little  sense  to  warrant  any 
statistical  treatment. 

However,  using  the  set  of  measured  n^/Nq  for  each  flame  and  knowledge  of 
Ttc,  an  effective,  averaged  value  of  the  Q/V  parameter  could  be  determined  for 
each  flame.  Using  these  three  values,  the  newly  determined  Tc  agreed  with  Ttc 
over  the  19  experimental  points  to  within  better  than  5Z.  This  means  an 
average  deviation  of  between  70  and  130K  depending  on  the  flame  temperature, 
which  varied  from  1400  to  2700K.  This  is  excellent  agreement. 

Why  does  the  calibration  method  work?  It  appears  that  the  effective 
value  of  Q/V  is  dependent  on  the  collisional  mix  but  not  much  on  the  tempera¬ 
ture.  It  Is  very  different  from  the  values  measured  directly  using  the 
downward  transfer.  Part  of  the  discrepancy  may  lie  in  the  dependence  of  the 
rotational  variation  of  the  various  collisional  rates  on  different  species. 
Also  recall  that  Q^/V  is  needed  for  the  analysis  of  N^/Nq  in  the  thermometry 
equation,  but  Q0/V  is  the  quantity  actually  measured  in  the  downward 
transfer.  From  low  pressure  measurements ,  ^  it  is  known  that  V  depends  on 
N'  in  v'  *  1  and  this  shows  up  in  flames  as  well.'*  Recent  measurements  in 
this  laboratory*^**^  show  that  Qc  is  N'-dependent.  The  dependence  is  true  for 
a  wide  range  of  collision  partners,  but  the  degree  of  variation  (always  a 
decrease  in  Q0  with  increasing  N')  Is  different  from  gas  to  gas.  No 
definitive  knowledge  exists  concerning  variation  of  Qj .  Even  though  the 
upward-transferred  rotational  distribution  is  near  thermal  (Section  III  B), 
perhaps  the  N'  dependence  of  these  quantities  over  differing  collisional  mixes 
is  simply  too  great  for  the  directly  measured  Q/V  to  be  applicable. 
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With  calibration  of  the  effective  value  of  Q/V  in  a  burnt  gas  collisional 
mixture,  very  good  temperature  determinations  have  been  made.  One  can  expect 
to  measure  temperatures  with  a  typical  precision  of  100K  in  an  exhaust  gas 
system  at  2000K,  and,  performed  this  way,  VETT  has  definite  promise  as  a  flame 
thermometer.  The  lack  of  a  fundamental  understanding  of  the  differences  bet¬ 
ween  the  values  of  (A/V)gff  and  the  measured  Q/V,  however,  suggests  extreme 
caution  in  setting  limits  of  applicability  of  VETT  to  various  flame  gas 
environments. 

Because  of  the  lack  of  correspondence  of  the  measured  Q/V  and  the 
(Q/V)eff  values  determined  by  upward-transfer  runs,  we  did  not  pursue  the 
proposed  three-channel  thermometry,  where  v'  ■  1  is  pumped  and  each  of  v'  -  0, 
1  and  2  is  observed.  This  was  in  spite  of  the  ability  to  discern  further 
upward  transfer  v'  ■  0  •*  1  -*•  2  (see  Figure  2).  Also,  in  view  of  the  lack  of 
fundamental  understanding  of  the  point  measurements,  imaged  thermometry  was 
not  attempted. 

There  are  other  possible  ways  to  measure  temperature  in  flames, 
o  n 

Cattolica  has  used  a  two-line  method,  in  which  the  same  upper  state  level  is 
excited  from  two  different  lower  rotational  or  vibrational  levels.  This 
avoids  any  problems  due  to  collisional  transfer  and  quenching,  but  does  so  at 
the  expense  of  accuracy  and  dynamic  range  in  the  temperatures  so  obtained.  If 
one  understood  the  relative  energy  transfer  and  quenching  rates  as  a  function 
of  N' ,  these  restrictions  could  be  relaxed  and  different  N'  levels  could  be 
pumped.  Much  of  the  pertinent  information  on  upper  state  OH  transfer  is 
contained  In  Refs.  5-7  and  15-19,  but  the  present  results  on  0/V  indicate  a 
greater  understanding  is  needed  in  order  to  use  such  an  approach.  Further 
research  along  such  lines  would  be  valuable,  because  a  two-line  method  pumping 
different  N'  could  be  implemented  with  our  recently  demonstrated  two-frame 
Imaging  method,  to  provide  an  instantaneous  map  of  the  temperature  field 
correlated  with  OH  concentration. 
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Figure  3.  Boltzmann  plots  of  relative  rotational  populations  in  the  v’-l  state 
following  upward  transfer.  Triangles  are  Fj  levels  and  circles  are  F2 
levels.  Top,  F^)  pumped  In  v'-O,  0.2  A  resolution,  T-2020K.  Middle 
Fj (5)  pumped,  0.8  A,  iwj9O0K.  Bottom,  F2(ll)  pumped,  0.8  A,  T-2130K. 


Figure  4.  Exemplary  thermometry  runs  as  stored  on  the  computer.  Plotted  Is 

Intensity  vs  wavelength  of  the  spectrometer,  which  is  operated  in  low 
resolution.  Top,  (0,0)  band,  initially  pumped;  bottom,  (1,0)  band 
populated  by  upward  transfer.  The  intensity  scale  for  the  lower  run  i 
*  factor  of  37.5  more  sensitive  than  in  the  upper  scan.  Integrated 
intensities  are  used  for  the  data  analysis. 


Figure  5.  Exemplary  Boltzmann  plots  of  rotational  populations  in  the  ground  state 
of  OH,  obtained  by  excitation  scans  and  used  to  determine  the 
rotational  temperature  TR.  Circles  and  triangles  denote  two  separate 
series  of  runs;  lines  represent  least-squares  fits.  The  corresponding 
temperatures  are:  circles,  1998K,  triangles,  2025K. 
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